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Abstract. We consider the graded S^-modules of higher diagonahy harmonic polyno- 
mials in three sets of variables (the trivariate case), and show that they have interesting 
ties with generalizations of the Tamari poset and parking functions. In particular we 
get several nice formulas for the associated Hilbert series and graded Frobenius charac- 
teristics. This also leads to entirely new combinatorial formulas. 



Our global aim here is to give an explicit description of the space IK^ of trivariate 
diagonal harmonic polynomials. We make interesting progress along these lines, 
and propose a combinatorial description of the graded Frobenius characteristic for this 
space. This description encompasses the "shuffle conjecture" of [15], and suggests new 
combinatorial identities linked to the combinatorics of the Tamari lattice. In fact it has 
lead us to introduce a new notion of r- Tamari poset, and these identities may be stated 
in terms of positive integer polynomials in the variable r. The first drafts of this paper 
has prompted other researchers to consider these combinatorial questions, and some have 
already been settled (see [5] and [6]). 

More explicitly, in the ring IZn '= C[X] of polynomials in the three sets of variables 
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"Kji is the subspace of polynomial zeros of the polarized power sums differential op- 
erators 

n 

P^{dX):=Y,dxpypzl (1) 

for a = (a, 6, c) running thought the set of vectors of norm \a\ := a + b + c, with 1 < 
\a\ < n. Here, dv denotes derivation with respect to v. It is clear from the definition 
that "Kn is closed under derivation. For one set of variables, say x = Xi, . . . "K^ is 
entirely described by classical theorems (see [19, section 3.6]) as the span of all partial 
derivatives of the Vandermonde determinant 

A„(x)=det o<,<n-i' 

It is of dimension n!, and isomorphic to the regular representation of for the action 
that permutes variables. The description of the case of two sets of variables took close to 
15 years to be finalized (see [17]), and much remains to be understood. It is of dimension 
{n + 1)^"-^. It has been observed by Haiman (see [16]) that the dimension of "K^ seems 
to be given by the formula 

dimJ{^ = 2^(r^ + l)^-^ (2) 
but, until very recently, almost no one had further studied the trivariate case. We have 
endeavoured to do so, considering more general spaces "Kn^ (see definition (12)) for which 
we have experimentally found out that 

dimlK^;) = (r + 1)^ (rn + 1)^-^ (3) 

Symmetric group action. We turn both IZn and "Kn into §^-module by considering 
the diagonal action of symmetric group on variables. Recall that, for a G Sat, the 
polynomial a • f{X) is obtained by replacing the variables in f{X) by 

2/(7(1) y(j{2) 

We denote by the monomial 

yA ._ ai ^ ^ ^ an . , , , bn . ^ci _ ^ Cn 

for a matrix of integers 

A = 

Writing Aj for the j^^-column of A^ we may define the degree of a monomial X^ to be 
the vector 




deg(X^) := 



in N^. Clearly the diagonal action preserves degree. The total degree tdeg(X^) of a 
monomial X^ is the sum of the components of deg(X^). An interesting subspace of "Kn 
is the space An af alternating polynomials in IK^. Recall that these are the polynomials 
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such that (J • f{X) = sign(a) f{X). It is also worth recaUing that, in the case of two sets 
of variables, this subspace has dimension equal to the Catalan number 

n + I \ n 

Another observation of Haiman (ioc. sit.) is that, in the trivariate case, we seem to have 

n{n + 1) \ n — 1 ) 

More generally, our experimental calculations suggest that the alternating component 
J^n of yi^n' has dimension 

d.m^!.-)^ (r + lW (r + l)^. + rY 

n(rn + 1) V ^ - 1 / 
It is intriguing^ that this formula bears a resemblance to a formula that appears in [?] . 

Degree grading. The space "K^ is graduated by degree, and its homogeneous compo- 
nents IK^^d (with d G N^) are S'^-invariant, hence is also graded. There corresponds 
direct sum decompositions 

— 3^n,d5 WiA J\ji — ^n,d- (6) 

As discussed in Section 2, the space "Kn is finite dimensional. In fact, we will see that 
the homogeneous component IK^^d is non vanishing only if 



(il + (^2 + ^3 < 



where d = (c/i, 6/2, ^3). We can thus consider the Hilbert series: 

(^) © 
^>^n(q) := dim(J{.,d) and yi,(q) := dim(yi,,d) (7) 

|d|=0 |d|=0 

of "Kn and An^ writing for c^c^c^ . The Hilbert series IK^(q) is a symmetric polyno- 
mial in gi, ^2 and ^3. In fact, it is always Schur positive as discussed in Section 3. This 
is to say that it expands as a positive integer coefficient linear combination of the Schur 
polynomials 5^(q) in the variables q. 

For example, we easily calculate that 3^2 affords the linear basis 

{1, X2 - xi, 2/2 - Z2 - zi}, (8) 

whose subset {x2 — ^i, 2/2 — ^2 — ^1} is clearly a basis of A2. We thus get the Hilbert 
series 

3^2(q) = 1 + (gi + ^2 + 93) = 1 + and yi2(q) = qi + q2 + qs = Si{q). 



■'^Thanks to the anonymous referee for this observation. 
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To calculate IK^(q) for larger n, we exploit the fact IK^ is closed under taking derivatives, 
as well as applying the operators (see Section 2) 

n 

1=1 

Here u, v stand for any two of the three sets of n variables: 

X = xi, . . . y = 2/1, ... ,2/^, and z = Zi, . . . 

It easy to verify (see Section 2) that the Vandermonde determinant A^(u) belongs to 
IK^, whether u be x, y of z. Using all this, we can calculate that 

J{3(q) = 1 + 2(^1 + ^2 + 93) 

+2 {ql + ql + ql) + ^ {qiq2 + qiqs + 92^3) 

+ql + ql + ql + qlq2 + qlq^ + qiql + qiql + qlq^ + 92^3 + 91^293 
by checking directly that we have the following respective bases Ba for each "Jin^d'- 
S300 = {A3(x)}, !B2oo = {a,,A3(x),c>,,A3(x)}, Sloo = {a2^A3(x),c>,,c>,,A3(x)}, 
So3o = {A3(y)}, S020 = {dyA^iy). dyAM). ^oio = {5,^, A3(y), dy,dyA^{y)}, 
Soo3 = {A3(z)}, Soo2 = R,A3(z),a,,A3(z)}, -Booi = A3(z),a,,a,,A3(z)}, 

= {^yxA3(x)}, -Biso = {£^yx^yxA3(x)}, '^,201 = {^zxAslx)}, 

Sl02 = {^zx^zxA3(x)}, S02I = {^.yA3(y)}, !Boi2 = {^zy^zyA3(y)}, 

Siio = {^yx5,,A3(x),Eyxa,,A3(x),4'x^A3(x)}, 

Sioi = {^.x5.,A3(x),E,,a,,A3(x),£;£)A3(x)}, 

Son = {^zy5,,A3(y),E.ya,,A3(y),Eiy^A3(y)}, 

Sooo = {1}, Sni = {^.x^yxA3(x)} 

Collecting previous observations, and doing some further explicit calculations, we get 
that 

^i(q) = 1, 
J{2(q) = l + si(q) 

J{3(q) = l + 2si(q) + 2s2(q) + Sn(q) + S3(q), (9) 
J{4(q) = 1 + 3si(q) + 5s2(q) + 3sii(q) + 6s3(q) + 5s2i(q) + siii(q) 
+5 S4 (q) + 4 S31 (q) + 3 S5 (q) + S41 (q) + S6 (q) • 

If we specialize one of the parameters (say ^3) to 0, we get back the graded Hilbert series 
of bivariate diagonal harmonics (of overall dimension (n + l)"""^) which has received a 
lot of attention in recent years (see [8, 10, 13, 14, 15, 18, 21]). In other words, the z- 
free (or x-free, or y-free) component of CK„ coincides with the "usual" space of diagonal 
harmonics (in the bivariate case, as previously considered in the literature). Hence, all 
of our formulas involving the parameters q specialize to known formulas by the simple 
device of setting one of the three parameters in q equal to 0. Evaluating IK^(q) at qi 
equal to 1, we clearly get the overall dimension of IK^. For this evaluation, the value 
of (9) agrees with formula (2). 
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Graded character. We refine the dimension analysis by taking into account the de- 
composition into irreducibles of the homogeneous components of IK^, see (6). This is aU 
encompassed into the graded Frobenius characteristic of IK: 

J{„(q;w) := J]qdj-^„^^(w). 

RecaU that the Frobenius characteristic J^v(w), of a §^-module V, is the symmetric 
function (in auxihary variables w = Wi^W2^ • • •) whose expansion in terms of the Schur 
functions Sx{w) records the multiplicity of irreducibles in V. This is to say that we have 

with the sum being over partitions of n, which classify irreducible representations of 
For reasons also discussed in [2], the expansion of IK^(q; w) in terms of the Sx{w) affords 
Schur positive coefficients in the 5^(q). This is to say that 

(q) 5a(w), nA,^ e N. (10) 

Ahn \ /i / 

Hence there are two kinds of Schur function playing a role here. To emphasize this, 
we denote those in the q- variables by a lower case "s", and drop the variables q. For 
example, we have 

J{i(w;q) = 5i(w) 

^2(w; q) 5'2(w) + 5iS'ii(w) 

3^3(w; q) = S^{w) + (^i + S2) 5^21 (w) + (^n + 53) S'iii(w) 

J{4(w; q) = 5'4(w) + (^i + ^2 + s^) S^i{v^) + (52 + ^21 + 54) 5'22(w) 

+ (^511 + 53 + 521 + 54 + 531 + 55) 5'21l(w) 
+ ('5111 + ^531 + 541 + Sq) S'iiii(w) 

3^5 (w; q) = S^{w) + (51 + 52 + 53 + 54) 5^41 (w) 

+ (^2 + 53 + S21 + 54 + 531 + 522 + 55 + 541 + 56) S'32(w) 

+ (511 + 53 + 521 + 54 + 2 531 + 2 55 + 541 

+ 532 + 56 + 551 + 57) 5^311 (w) 

+ (521 + 54 + 531 + 522 + 5211 + 55 + 2 541 + 532 

+ 5311 + 56 + 2 551 + 542 + 57 + 561 + 53) 5'22l(w) 

+ (5111 + 531 + 5211 + 2 541 + 532 + 5311 +56 + 2 551 

+ 542 + 5411 +533 + 57 + 2 561 + 552 + 58 + 571 + 59) 5^2111 (w) 

+ (5311 + 542 + 5411 + 561 + 5511 + 543 

+ 571 + 562 + 581 + 5io) S'iiiii(w) 

(11) 

In these expressions, the Schur functions 5^ = 5^(q) are to be evaluated in the three 
parameters q = qi^q2^Q3- The Schur function aSa(w) act as placeholder for irreducible 
representations of S^. Thus the multiplicities of the irreducible representations of 
in IK^, classified by the partitions /i of n, appear as the coefficients of ^^^(w) through 
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the evaluation of (11) at gi = 1, ^2 = 1 and qs = 1. In particular, the value obtained 
this way, as coefficient of 5'ii...i(w), agree formula (4). Furthermore, the expressions in 
(9) may be directly calculated from (11) by replacing each S'a(w) by the number of 
standard tableaux of shape /x (see [1]). 



Generalized harmonics. Using notions further discussed in Section 2, all of these 
considerations generalize to the spaces of higher diagonal harmonics 

!K« := n (^;-%)^ (12) 

where An (resp. X^) stands for the ideal generated by alternating polynomials (resp. 
invariant) polynomials in 7^, and the orthogonal complement is with respect to the 
§^-invariant scalar product defined by 

iX-,X-):^h' 'f;^ = ^' (13) 

I otherwise. 

Here A\ stands for the product of factorials: 

A\ := ai! • • - a^! • 6i! • • -bj • ci! • --cj. 

(r) 

When r is even, the diagonal action on 'Kn is twisted by the sign representation. For 
n = 2, the space 3^2^^ decomposes as a direct sum of its isotypic invariant component, 

(r) (r) 

^2 , and alternant component, A2 : 

:K^^=3^^®aP, (14) 

which afford the respective bases 

Sa^^^ = {{x2-xiy{y2-yi)\z2-ziY\a + b + c = r-l}, and 
S4^^ = {{x2-xir{y2-yi)\z2-ziY\a + b + c = r}. 

It follows that 

3^?^(w;q) = Sr-lS2{^^r) + SrSu{w), 
from which we get the Hilbert series 

^2^^(^) = '^r-l + J2^^(q) = ^r-l, and ^2^^ (q) = Sr. 

In particular, 

dimJCr^ = (r + lf, dnn3i^ = '±±^, and dimA^^ = ^^^^±^^^^. 
which agree with formulas (3) and (5). 

The above formulas, which clearly generalize (2) and (4), are both encompassed in the 
following generalization of a Conjecture of [15], in conjunction with a second conjecture 
below. 
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Conjecture 1. The graded Frobenius characteristic ofKn"^ is given by the formula 

J{W(w;gi,g2,l)= E E ^^(^i)92''^^^^Qco(/)(w), (15) 

where (5 run over the set of r-Dyck paths, on which we consider the r-Tamari order 
(see Section 5)^ andVJ-'{/3) denotes the set of ^^parking functions of ^^shape^^ (5. 

To finish parsing the right-hand side of (15), recall from [15] that (5co(/)(w) stands for the 
fundamental quasi-symmetric polynomials indexed by the composition co(/) associated 
to a parking function / as described in Section 5. We also recall in Section 5 other 
relevant notions pertaining to Dyck paths, parking functions and the Tamari poset, 
including the definition of the "dinv" statistic. For a given r-Dyck path /3, let us denote 
by i^p (q) the polynomial 

a</3 

that g-enumerates the r-Dyck paths a that are smaller than /3 in the r-Tamari order. 
Such paths are weighted by where d{a^ /3) stands for the of length of the longest 

chain going from a to ;5 in Vn"^ . We write i^^^ for i^^^(l). Since rf(a, /5) = if and only if 
Of = /3, we have ^^^^(0) = 1. Thus Conjecture 1 agrees with the shuffle conjecture of [15] 
at gi = 0. It also follows from results of [ ] that the specialization at gi = ^2 = 1 of 
formula (15) gives 

J{M(w; 1,1,1)= V 4^)eeo(/3). (16) 



Here e^^...^^ := 6^ • • • ^rfc? where stands for the usual r^^ elementary symmetric function 
(see [22]). It may be worth underlying that the right hand side of this last equation is 
simply, up to twisting by the sign, the Frobenius characteristic of the action of the 
symmetric group on the set of pairs (/, Qf), where / is a r-parking function and a is a 
r-Dyck path lying below the shape of /. See section 5 for more on this. Calculations 
suggest that we have the following explicit formula for this Frobenius characteristic. 

Conjecture 2. 

JC«(w; 1, 1, 1) = (-l)"-^W(rn + 1)^^- H ^) f ^^^^^^ ^^^^ 

Recall that, for a partition A having di parts of size it is usual to denote by zx the 
integer 

Zx := l^Mi! 2^2^21 ••• n^^rfn!. 

It is interesting to note that the special case r = 1 of this last formula has also been 
conjectured independtely by Loktev (see [21]). It follows from (17) that we have an 
explicit formula for the enumeration of pairs (/, o^), as above, that are fixed points for 
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the action of a permutation having given cycle type. In generating function format, this 
may be stated as 

E E ^^^co(,) = -^exp('E(- + l^ (18) 

n>o^^^(r) V ; V / / 

where the e-basis has been turned into the /i-basis by the removal of signs. 

This also implies that the dimension of IK^^ may be described in two different ways, 
giving rise to the following elegant formula 



A first draft of this paper prompted a collaborative eflFort between Bousquet-Melou, 
Chapuy, and the second author, to give a direct proof of this new combinatorial identity 
(see [6]). 

On another note, we will recast formula (17) using the following calculations. We start 
with establishing the formal power series identity 

\ k>l ^ ^ / n>l 

where we use the Pochhammer symbol notation 

(^)(^) :=u{u + l){u + 2) (u + m-l). 

To check that the above equality hold, we may proceed as follows. Let us denote by 

Z = Z{t) the series corresponding to the right-hand side of (20), at a = 1. Then 
Equation (20) is equivalent to 

E n - f . pi) 



k>l 



On the other hand, for 6 G N, we see directly that the series Z is algebraic and satisfies 
the equation 

Z^-^ = {Z- t)\ (22) 
with initial condition 2^(0) = 1. But (22) implies that we have 

^ - ^ (23) 

and identity (21) follows from this. 

We now use Formula (20), setting a = r n -\- 1 and 6 = r + 1, in conjunction with the 
symmetric function identity 



exp 



^p,(u)^) =^e.(u)r, (24) 



k 

k>l / n>0 



TRIVARIATE DIAGONAL HARMONICS 9 

to "evaluate" the elementary symmetric function e^(u) for the plethystic substitution 
(see [1]) 

p,(u)^(rn + l)(^^'^Y^^)- (25) 

Recall that this means that we should expand en{u) in terms of the power-sums p/c(u), 
and then apply (25). In view of (20), we get that 

This calculation allows us to reformulate (16) in terms of the monomial symmetric func- 
tions mx{w) by applying the substitution (25) to the Cauchy identity (see [22, page 
65]): 

J](-l)"-^WpA(u) -Pa(w) = J2 ^aH m,(w). (27) 



We thus get 



5](-l)-^W(rn + l)^Wn Y^^) fpA(w) 



E/ .\f(\)TT 1 fr((r + l)n - k + 1)\ , , 



(28) 



simply by applying the substitution (25) in the left-hand side of (27), and (26) in the 
right-hand side. 

Formula (17) can thus be rewritten in the form 
Mi:«(w; 1.1,1) = + n f'"'' + - * + m.(w), (29) 

Ahn keX ^ ^ 

In particular, the multiplicity of the alternating representation in J-Cn is the coefficient 
of mii...i(w) in this last expression. In view of (16), this suggested that we have the 
following generalization of a formula of Chapoton (see [7], in the case r = 1) for the 
number of intervals in the r-Tamari poset: 

^ ~ n{rn + l)[ n-1 J' 

A few months before this writing, the first author conjectured that this last equality 
should hold. This new combinatorial identity has since been shown to be true in [5]. 

In view of Proposition 3, the above formulas and conjectures imply that the following 
combinatorial identity should also hold 

4^^X(co(/5) = ll...l) = ^"'^ (31) 
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where X{P) is equal to 1 if P is true, and otherwise. A direct bijective combinatorial 
proof of this identity seems to follow from the approach in [o]. Moreover, it appears that 
(18) may be amenable to a Polya theory extension of the approach of [ ] and [5]. In 
particular, this will reduce Conjecture 2 to Conjecture 1, as well as giving a joint proof 
of (19) and (30). 

2. Trivariate diagonal invariants and alternants 

Invariants and scalar product. For our purpose, some interesting properties of trivari- 
ate diagonal invariants and alternants need to be recalled. Much of what is discussed 
here is a straightforward generalization to the context of three sets of variables of results 
of [12] for the bivariate case. They are included here to help make our presentation 
understandable on its own. Omitted proofs are entirely similar to the bivariate case. 
A particular case of a result of Weyl, in [24], establishes that the subring of trivariate 
diagonal invariants is generated by the polarized power sums 

n 

P.(X):=5^X;, (32) 

for a = (a, 6, c) running thought the set of all N- vectors such that 1 < |o^| < n. 

To make better sense out of Definition (12), we need to discuss a few notions concerning 
trivariate diagonal alternants. Let us first observe that definition (13) is equivalent to 

{fiX),g{X)) = fidX)giX)\^^^, (33) 
where p{X)\^_^ := p{0) is the constant term of p{X). 

It is a well known fact that, with respect to this scalar product, the orthogonal com- 
plement of an ideal / = (/i, . . . , /at), is simply the set of polynomial zeros of the 
differential operators fi{dX). Indeed, let g{X) be in and consider the leading term 
(for any suitable term order) of 

MdX)giX) = cX^ + ..., 

if any. Since X^ fi{X) lies in /, we must have {X^ fi{X)^ g{X)) = 0, but this means 
precisely that c = 0. Hence we must have fi{dX)G{X) = 0. In particular, Weyl's result 
implies that 

To show that the spaces we consider are finite dimensional, the following is very useful. 
Proposition 1. All monomials X^^ of total degree larger than (2) lie inX^. 

Alternants. A set of linear generators for trivariate alternants is easily obtained by 
applying the antisymetrization operator 

i?±(X^):= J]sign(a)a•X^ (34) 
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to any monomial X^, with ^4 a 3 x n matrix of non-negative integers. Evidently the 
product of an alternant by an invariant is yet again an alternant. We get a linear 
basis Aa{X) := of the module of alternants over the ring of invariants, 

by choosing matrices A having all their columns different and appearing in decreasing 
lexicographic order (from left to right). Any symmetric operator in the ring C[X^dX] 
sends alternants to alternants. In particular, this is the case for the operators £'uv, 
and Pa{dX). A degree argument shows that Pa{dX)AA{X) = 0, whenever A is an 
order-ideal under component comparison of columns: 

B <C, iff {b, < q). 

For instance, this is the case for the Vandermonde determinants. Moreover, a simple 
direct calculation of the commutator of the operators E'yx and Pa{dX)^ with a = (a, 6, c), 
gives that 

[PaidX),E^''2] = aPp{dX), (35) 
where /3 = + (fc, — 1, 0). It is evident that similar results holds for other possible choices 

(k) 

of variable sets in Euv • Putting all this together, we deduce that 
Proposition 2. If A e N^><^ is an order-ideal, then 



{dX^)n{AA{X)) G (36) 



for any A, and any composition Q of the operators Eu^ - 



The "Operator Theorem" of [18] states that Proposition 2 entirely characterizes the 
space "Kn in the bivariate case. More precisely, we have the following. 

Theorem (Haiman). The space of diagonal harmonics for the bivariate case is the small- 
est vector space containing the Vandermonde determinant A^(x)^ which is closed under 

(k) 

taking partial derivatives and applications of the operators Ey^ . 



For the trivariate case, experiments suggest that the analogous statement should hold. 
However, the methods employed by Haiman to settle the bivariate case do not seem to 
generalize to the trivariate case. 

Another result along these lines is that, for the bivariate case, we can calculate the entire 
space aP (resp. J^'-)) by successive applications of the operators £'yx , starting with the 
r*^-power (resp. (r — l)*^-power) of the Vandermonde determinant A^(x). However, we 
do not know how to explicitly construct a basis in this manner. All, this also appears to 
hold in the trivariate case. In particular, we have that 

Proposition 3. For all n, and all r >\, the spaces and An~^^ coincide. 
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3. Action of the general linear group 

To make apparent a special feature of both the Hilbert series and the (graded) Frobe- 

(r) 

nius characteristic of the spaces 'Kn , we need to consider it as a polynomial GLa-sub- 
representation of TZn = C[X], for the action 

f{X)^ f{MX), for MeGLs. (37) 

Clearly this action commutes with the action of Sn- Recall that the character of a GLs 
representation is the symmetric function of the parameters q = qi^q2^Q3 obtained by 
calculating the trace of the linear transformation 

Q*{f{X)) :=f{QX), 

where Q is the diagonal matrix Q = [qi^q2^Q3]- Observe that a polynomial f{X) is 
homogeneous of degree d if and only if 

/(gX) = q''/(X). 

Since an homogenous subspace V, of 7^^, is simply a subspace that affords a basis B of 
homogeneous polynomials, it is Ga-invariant and its character V(q) is 

V(q) = 5]q^^^(^). 

feB 

This is precisely what we have called previously the Hilbert series of V. It is well 
known that the characters, of irreducible polynomial GLa-representations, are the Schur 
polynomials (evaluated at gi, ^2, Qs)- Hence, V(q) expands as a positive integer coefficient 
linear combination of the 5^(q), since the coefficients correspond to multiplicities of 
irreducibles. 

Since the two actions commute, we get a double decomposition into irreducibles, for both 
E>n and GI/3, of any subspace that is stable under the two actions. This is the reason 
why the graded Frobenius characteristic of such a space expands as a positive sum of 
products Si^{q)Sx{w). 

4. The w^hole story for the case n = 3 

Using observations and results of [2], we can calculate explicitly J{3^^(w;q) in full gen- 
erality. The crucial observation is that the only Schur function 5^(q) that can occur in 

the expansion of :K^i\w ; q) are those indexed by partitions having at most two parts. 
We can thus deduce the trivariate case^ expansion from the known bivariate case. 

Recall that, in the bivariate case we have J{f^(w;g,t) = V(5iii(w)), where V IS an op- 
erator characterized as follows. As usual (see [22]), we denote by A ^ /x the dominance 
order on partitions, and /x' stands for the conjugate partition of //. The integral form 



'In fact, this particular calculation holds for any number of set of variables. 
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Macdonald polynomials i?^(w; t), with fi a partition of n form a linear basis of the 
ring A, of symmetric functions in the variables w, and are characterized by the equations 

(i) {S>,{w),H,[{l-q)w-q,t])=0, if A ^ ^, 

(ii) {Sxiw),H^[il-t)w;q,t])=0, if A ^ and (38) 

(iii) {Sniw),H^{w;q,t)) = 1, 

involving the usual "Hall" scalar product on symmetric functions (for which the Schur 
functions are orthonormal) , and some plethystic notation (see [!]). The linear operator V 
is characterized by the fact that it affords the Macdonald polynomials as eigenfunctions: 

ViH,iw;q,t)) = H,iw;q,t) J] q'^t\ 

This presentation may seem a bit terse, but for our purpose we only need the case n = 3 
for which an explicit description of V is given below. Using this description, we will infer 

(r) 

the graded Frobenius of "K^ for the trivariate case, out of the simple knowledge of the 
bivariate case. To this end, we consider the linear operator corresponding to the matrix 

/O 522 ^32 \ 

V= -521 -531 , (39) 

\1 51+52 511+53/ 

with respect to the basis {aS3(w), 6^21 (w), aSiii(w)}. Then^, 
Proposition 4. For the trivariate case, we have the explicit formula 

^r(w;q) = V'^(5ni(w)). (40) 

We may calculate from (39) that 

Lemma 1. The characteristic polynomial of the operator V , over the ring A(^)(q)^ is 

- (53 - 521 + 5ii)V2 + (541 + 533 - 532)V - 544. (41) 

One may observe that 5ii(q) = ^1^2 + ?i?3 + ?2?3 is one of the eigenvalues of V. The 
associated eigenspace is spanned by the symmetric function 

5'3(w) + {qi + (72 + ga) 5^21 (w) + ((71(72 + (7i<73 + <72<73) 5'iii(w), 

which is a three parameter version of a Macdonald polynomial. Wether one can develop 
a theory of these polynomials for all n remains an open question. One would expect 
these to appear as eigenfunctions for a suitable generalization of V. 



^All calculations are to be done over the ring A^^^(q) of symmetric functions modulo the ideal gen- 
erated by the for which ji has at least three parts. 
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In view Proposition 4, Lemma 1 leads to a recursive approach for the calculation of the 

(r) 

graded Frobenius characteristic of any of the spaces IKg , based on the initial values 
J{f(w;q) = 5ni(w), 

^3'^(w;q) = SsM + {s2 + Si)S2iM + {ss + Sn)Sin{w), 
^^(w; q) = (511 + Ss)Ss{w) + (521 + S4 + ssi + 55)^21 (w) 

+ (^52,2 + ^541 + 56)S'iii(w). 

Calculations are to be done in the ring A(^^(q) of symmetric functions spanned by Schur 
functions indexed by partitions that have at most two parts. In order to describe the 
resulting formulas, let us introduce the notation 

Tnk '•= + 5^-2,1 + . . . + Sn-2k,k' (42) 

After calculations, we get the following. 

Proposition 5. The space IKg^^ affords the explicit graded Frobenius character formula 

3^^^(w;q) = T3(^_i),^_iS'3(w) + (r3^_2,r-i +r3^_i^^_i)5'2i(w) + r3^,^S'iii(w). (43) 

In particular^ the Hilbert series of the alternating part A^^"^ of this space appears as the 
coefficient of Siu{w) in this formula^ i.e.: 

^f^(q) = S3r + S3r-2,1 + • • • + 5^,r. (44) 

For any a > 6 > 0, the Schur function Sab = Sab{Qi^ 92,93) specializes to 

^(6+l)(a + 2)(a-6+l), 

when one sets 9i = 92 = 93 = 1- From this, it follows that 

T,fc(l, 1, 1) = ^ (fc + 1) {k + 2) - (13 + lOn) A; + 3 (n + 1) (n + 2)) 

We may then check directly that formula (43) specializes at gi = ^2 = 93 = 1 to the case 
n = 3of (17). 

5. Dyck paths, Tamari posets, and parking functions 

Generalized Dyck paths and Tamari posets. Recall that a r-Dyck path of height 
n is a path in N x N, that starts at (0,0), ends at (rn,n), and stays above the line 
X = ry. For instance, we have the 2-Dyck path of Figure 1. Such a path may be entirely 
described by a sequence a = (ai,a2, . . . ,a^), where (a^,i) corresponds to the leftmost 
point lying on the path at height i. I other words, is the number of "boxes" lying to 
the left of the path at height i. In this manner, a is a r-Dyck path if and only if 

(1) < ai < a2 < . . . < a^, 

(2) for each i, we have < r{i — 1). 
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(10,5) 



(0,0) 

Figure 1. The 2-Dyck path encodes as 00367. 

(r) 

We denote by Vn the set of r-Dyck paths of height n. For example, 

^ = {000, 001, 002, 003, 004, Oil, 012, 013, 014, 022, 023, 024} 

It is weU known (see [9, 11]) that the number of r-Dyck paths of height n is given by the 
Fuss-Catalan number: 

^ rn + 1\ n J 

We say that (a^, a^+i, . . . , ak) is a primitive subsequence, of a r-Dyck path (ai, . . . , a^), 
if 

(1) aj — ai < r{j — i) for each i < j < and 

(2) either = n, or a/^+i — > r{k + 1 — i). 

For each i, there is a unique such primitive subsequence. It corresponds to the portion 
of the r-Dyck path that starts at (a^, i) and ends at the "first return" of the path to the 
line of slope 1/r passing through the point (a^, i). Whenever i is such that a^_i < a^, we 
set a < /5, where 

/3 := (ai, . . . , a^_i, — 1, . . . , a/^ — 1, a^+i, . . . , a^), 

with {ai^...^ak) equal to the primitive subsequence starting at a^. The r-Tamari 
poset"^, on the underlying set Vn\ is the order obtained as the refiexive and transi- 
tive closure of this covering relation a < Its largest element is the path 00 • • • 0, and 
its smallest is the path for which = r {i — 1). 

For a given r-Dyck path /3, we have already considered the number i^^^ of r-paths that 
smaller than /5 in Vn"^: 

t^P := #{« G I « < /?}, 

and denoted by /5) the length of the longest chain going from a to /5 in Vn\ Observe 
that, when a is the smallest element of Vn\ d{a^/3) coincides with the usual "area" 

^As far as we know this poset has not been considered before. 
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Figure 2. The 1-Tamari poset for n = 4. 



a r-Dyck path /5 we associate the composition co(/3) whose parts are the length of runs 
of 1 in /5. Thus co(00112) = 221. 



Parking functions. For a sequence of positive integers / = (/i, /2, • ■ ■ , /n) (we also 
write / = /1/2 ■■■/«), let 

f "(/)\ _ /«1 «2 • • • C 



(45) 



be the lex-increasing rearrangement of 

1 2 ... n 

Jl f2 ■ ■ ■ fn^ 

ordering first with respect to the bottom row. This is to say that bi < fej+i, for all 
1 < i < n, and that Oj < a^+i when hi = bi+i. One says that / is a r-parking function 
if /?(/) is such that bk < r{k — 1), for all k. Perforce, /3{f) is thus a r-Dyck path. It is 
said to be the shape of /. 

For example, with r = 2, we have the following 49 r-parking functions of length 3: 

000 001 002 003 004 010 Oil 

012 013 014 020 021 022 023 

024 030 031 032 040 041 042 

100 101 102 103 104 110 120 

130 140 200 201 202 203 204 

210 220 230 240 300 301 302 

310 320 400 401 402 410 420. 
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The set of r-parking functions of length n, denoted by VJ^^^\n)^ has cardinahty (rn+l)^~-^ 
(see [25]). For a r-Dyck path /5, we further denote by VJ-'{/3) the set of all parking 
functions of shape /3. Observe that these are automatically r-parking functions. Clearly 
we get a permutation action of symmetric group on each of the VJ-'{/3) (by permutation 
of values). The corresponding Frobenius characteristic is easily seen to be equal to 
^co(/3)(w), where hc^...cf^{w) := /ici(w) • • • hcj^{w) is the complete homogeneous symmetric 
function associated to the composition co(/3) = Ci - • -Ck- In particular, the number of 
parking functions of shape /3 is given by the multinomial coefficient 

#V:F{P)=( r. (46) 

Vco(/3)y Ci\C2\"'Ck\ 

Formula (46) reflects the fact that any representation of S^, having /ic(w) as its Frobenius 
characteristic, must be of dimension (^). From this we get the identity 



(rn + 1)^-^ = 




Recall that it has been shown in [18] that the (ungraded) Frobenius characteristic of the 
z-free component of IK^ is equal to the character of the S^-module VJ^^^\n)^ twisted by 
the sign representation. From the point of view of Frobenius characteristics, this twist 
turns hc{'w) into ec(w). Hence, one consequence of the results in [ ] may be stated as 

J{M(w; 1,1,0)= eco(/3)(w). (48) 

The composition co(/) = C1C2 • • • c/^, associated to a r-parking function /, is simply the 
composition that encodes the descents of the permutation a{f). 

For the purpose of stating Conjecture 1, and discussing potential extensions (see next 
section), let us recall (from [ ]) how to calculate the "dinv" statistic of a parking 
function. Given a r-parking function /, let and bi be as in (45), and let Ci := ri — hi. 
For 1 < i < j < n and < < r — 1, we say that a triple (i, d) is a L^-inversion of / 
if either 

• Ci — Cj + d = and < aj , or 

• 1 < Q — Cj + < r — 1, or 

• Ci — Cj + d = r and > aj . 

Then dinv(/) is simply the number of D-inversions of /. Notice that this actually 
depends on the value of r. 

6. Final remarks 

A missing component in (15) is a "statistic" iy{f, a) that would account for the behavior 
of the third parameter ^3, with respect to pairs (/, a), for / G VJ-'{(3) and a < (3 in Vn\ 
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Such a statistic would give a complete combinatorial description of JC^T (w; Q'l, 0^2, ^3) in 
the form 

:KM(w;gi,g2,g3)= E E T.^'^^^'^^''^^^^^ (49) 

We would expect u to be such that a) = if and only if a is the smallest element 
of Vn\ Then, at qs = 0, formula (49) would specialize precisely to the Conjecture of 
Haiman et al in [15]. However, such a statistic still has to be found, even in the case 
r = 1. 

Much of what is discussed here seems to hold for generalized permutation groups G(r, n). 
This will be explored in upcoming work. 

Thanks. We thank Mark Haiman for insightful suggestions to the effect that we should 

(r) 

extend our approach to the case of the r-spaces 
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